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'\0 '. Abstract 

o 

^D ' In this paper we revisit thermodynamic constraints on primordial black hole (PBH) formation 

{^' in the early universe. Under the assumption that PBH mass is equal to the cosmological horizon 

^: 

^ , mass, one can use the 2nd Law of Thermodynamics to put a lower limit on the PBH mass. In 

^O ' models of PBH formation, however, PBHs are created at some fraction of the horizon mass. We 



> 

o 
o 

o 



1 



show that this thermodynamic constraint still holds for sub-horizon PBH formation. 
PACS numbers: 98.80.Cq, 04.70.Dy, 97.60.Lf 



I. INTRODUCTION 

The primary utility of primordial black holes (PBHs) [l|, |2| in cosmology today is their 
absence; their lack of observation allows one to put a limit on the size of density perturbations 
in the early universe J3|,y. Light {Mbh < lO^^g) PBHs will evaporate by the present day 



due to Hawking radiation emission 



have observable consequences 



. The injection of that radiation into the universe can 
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unless the PBH number density 



is sufficiently small. The PBH lifetime is related to the PBH mass, and the PBH mass is 
proportional to the horizon mass at which it formed. We therefore can use the lack of PBH 
emission at a certain time (the time of Big Bang nucleosynthesis until the present day) to 
constrain PBH formation at a much earlier time (from matter-radiation equality back to the 
time of reheating). 

Of the steps involved in this, the most uncertain is the relation between the PBH mass 
{Mbh) and the horizon mass {Mh) at which they formed. The mass of a primordial black 
hole formed at a time t is 

Mbh = JMh, (1) 

where the horizon mass Mh is defined through 

Mh = ^7rpR% (2) 

with the horizon size Rh = H~^ and the Hubble parameter H given by the Friedmann 
equation 

H- = ^Gp. (3) 

The function / is the fraction of the horizon mass that goes into the PBH. An analytic 
calculation in Carr [14] gave / = w^^^, where w = p/ p is the background equation of state. 
For a radiation dominated universe, w = 1/3 and so / ~ 0.2, close enough so that most 
authors simply assumed / ~ 1. The upper bound on / for a radiation dominated universe 
is unity, though it is lower for backgrounds with stiffer equations of state |l5|. 

The formation of a black hole changes the entropy of the universe. In the early uni- 
verse, the formation of PBHs leads to isocurvature (entropy) as well as adiabatic (density) 
aturbations |lq. The thermodynamic constraint on PBH mass, first considered in Lee 
, derives from the 2nd Law of Thermodynamics: that the entropy of a closed system 
can only increase. There, the horizon-sized region of the universe where PBH formation is 



happening is the closed system (in that it is initially causally disconnected from the rest of 
the universe). The entropy change we consider is that between the initial state radiation 
and the final state PBH-radiation system. For PBH creation exactly at the horizon mass 
(/ = !), there is no final state radiation and it is straightforward to compute the difference 
in entropy. The initial state consists of a perturbed FRW universe filled with a perfect 
radiation fluid of known density (and thus known horizon mass) . The final state consists of 
a PBH equal to the horizon mass with temperature and entropy given by 

K, 1 

^^^ " 2^ " SttGMbh' ^^^ 

Sbh = ^ = ^^GMIh, (5) 

where A and k are the PBH surface area and surface gravity, respectively. That the PBH 
must have more entropy than the initial radiation can be turned into the limit Mpbh ^ 
0.37Mp |l2| [2l|- We go through the derivation of this limit in a later section. Note, however, 
that the Lee bound is coincident with the implicit quantum "bound" at the Planck mass 
where quantum gravity becomes important. There we are no longer certain about the physics 
of black hole creation. For a horizon sized PBH, the thermodynamic bound is not higher 
than this quantum bound. It was not known is this continues to be true for sub-horizon 
PBH formation, which is the more physical situation. We prove in this paper that it does. 

II. SUB-HORIZON PBH FORMATION 



It was observed in 



ISl ligj l that gravitational collapse near threshold exhibits self-similar 



scaling behavior. The consequence of this for PBHs is that / is not a constant as shown 
above but rather depends on the size of the perturbation. For a PBH formed from an 
overdensity 6 = {p — p)/p being above a threshold 6c, the fraction is given by the power law 

f = K{6- 5c)\ (6) 

where k ?« 3, 7 ~ 0.37 and 5c ~ 2/3 for PBH formation during the radiation dominated 
era. Formally / vanishes as 5 ^ 5c- 



Numerical simulations of gravitational collapse [20|, |21|, |22| verify this power-law behavior 



though the values of the parameters (7, k, (5c) differ shghtly for different shapes of the initial 



TABLE I: Collapse parameters for different density perturbation profiles. 



Reference 



Gaussian 



Mexican hat 



Polynomial 



Niemeyer and Jedamzik [20] 



7 = 0.34 
5c = 0.70 



7 = 0.36 
5c = 0.67 



7 = 0.37 
5c = 0.71 





K = li.y 


K = 2.9 


K = 2.4 


Hawke and Stewart [211 


7 = 0.35 


- 


- 




5c = 0.4 


- 


- 


Musco et al. [22] 




7w 0.36-0.37 






5c = 0.71 


5c = 0.67 


5c = 0.71 



density perturbation. The profiles most common in numerical studies are a Gaussian, a 
Mexican hat and a sixth order polynomial (see Figure 1 of J20| for definitions and shapes), 
and we list the parameters in Table HI While the values of k and 5c vary by more than 
a factor of 2, the slope 7 is within 1% of the value 7 = 0.3558... computed analytically in 
Koike et al. [2^. 

The results of J2l| further suggest that / reaches a minimum value at ~ 10~^'^. Exactly 
how and why this occurs is not yet known, as this result has not yet been seen in other 
numerical simulations. If further studies show this bound to be a numerical artifact, then it 
remains that a PBH may have an arbitrarily small mass compared to the horizon at which 
it formed. We show in this short paper that there is indeed a lower bound on / coming from 
entropy constraints. 



III. ENTROPY CONSTRAINT 



For subhorizon PBH formation (/ < 1), one must take into account the entropy of the 
radiation outside of the event horizon but within the cosmological horizon. Using entropy to 
constrain PBH formation then becomes more subtle as the PBH is now no longer an isolated 
region. The entropy of the final state radiation is initial condition dependent {e.g., the shape 
of the perturbation as discussed in the previous section) and will have contributions from 
shocks. 

To be as general as possible, we do not attempt to account for the total amount of entropy 



in the horizon volume during collapse. Instead, we consider the region where a PBH is formed 
at the exact moment of PBH creation. During the infinitesimal time spanning the moment 
where a high radiation density region formed an event horizon and thus a black hole, there 
is a finite change in the total entropy. Any change in entropy due to an outgoing finite 
radiation flux, on the other hand, will be infinitesimally small. 

Our bound comes from comparing the entropy of the region immediately before and after 
PBH creation. We begin with the final state, where the PBH has a mass and entropy as 
given in Equation (jSJ. Our region of interest is then the sphere of radius Rbh = '^GMbh 
and volume Vbh = ^t^R'bhI'^- We assume throughout that the radiation is a perfect fluid, 
so that its density profile completely determines its behavior. We can define an average 
density of the radiation before PBH formation, 

^ Mbh_ ^ _3_Ml_ ^ p_ 
P* Vbh 2l7rM|^ P' ^> 

where the last equality is obtained from using Equations (P-El). 

The entropy density for a perfect radiation fluid of density p, pressure p and temperature 

T is given by 

Sraa='-^ = -^ = -('-^\'\''\ (8) 

T 3T 3 V 30 y '^ ^ ' 

where we have used the density temperature relation for a perfect radiation fluid with g^ 

effective relativistic degrees of freedom J3l| , 

P=^9.T^ (9) 

The total entropy inside the region of volume Vbh is approximated by 

Sraa=U^-4-\'' pTVbH. (10) 



3 V 30 
The condition that Sbh > Srad is 



^^ > _ /il^~o.39 / ^* (11) 



Mp - 27 V Stt ■ V 106.75' 
Or equivalently a bound on /: 



f>—9j—] ■ (12) 

•' - 405^ \Mp ' ^ ' 



This result verifes that of Lee [17| on sub-horizon scales: that the Planck mass remains 
the thermodynamic lower bound on PBH mass. Pushing this bound higher requires changing 



(7*, which is determined by the model of particle physics one uses. Including only Standard 



Model particles, g^ = 106.75 at T > 300 GeV 



24| . This will roughly double in SUSY models, 



ihough since the bound oc ^/g^, this has only a small effect. Only in a Hagedorn-type model 

2^ where the number of states increases exponentially with temperature could the bound 

increase appreciably. However, in that case the background cosmology we have assumed is 

likely to be invalid. The same is true for models of extra dimensions, where the fundamental 

scale of gravity is different from Mp. The formation and evaporation of PBHs is different 

26[ in that case and we do not consider it further here. 

IV. CONCLUSIONS 

We note that this bound might not be saturated in any astrophysical situation of PBH 
creation. This is because we have assumed throughout that PBHs form only at a single 
time (and thus single horizon mass). As PBHs will form for a range of horizon masses, the 
mass at which a PBH is formed (at a particular point in space) is given by maximizing the 
expression for Mbh over both S and M^- Using an excursion set formalism. Green and 
Liddle [27| showed that PBHs still form at order the horizon mass for the both the case of a 
blue perturbation spectrum and a spike at a specific scale (assuming a number density that 
saturates the current constraint). It is possible, however, that the result breaks down (so that 
we have / <^ 1) for a different type of spectrum and/or for a smaller number density of PBHs 
produced. However, in the limit that the PBH number density becomes astrophysically 
interesting (approaches the current upper bounds), there is still a tight correlation between 
PBH mass and horizon mass J28|. Our bound, therefore, would be satisfied by a large margin 
for any theoretically observable PBH population as a whole, but could be by a smaller margin 
for any individual PBH. 
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APPENDIX A: SELF-SIMILAR GRAVITATIONAL COLLAPSE AND CURVED 
SPACE EFFECTS 

To simplify the analysis in this paper we have assumed flat (Minkowski) space and uni- 
form radiation density, neither of which are strictly accurate in PBH formation. It might, 
however, be possible to exploit the self-similar nature of the collapse to obtain a more ac- 
curate calcuation of the initial state radiation. Self-similar perfect fluid solutions have been 
extensively studied [29|. The entropy calcuation including both curved space and using 
a proper density profile would be straight-forward though tedious. A self-similar solution 
asymptotic to the FRW universe (which is our initial conditioii) maintains a fiat density 



profile at small radii, which is verified in numerical simulations 211. The drawback to this 

n 

method, as first noted by Evans and Coleman [13] , is that the collapsing fluid only maintains 
the self-similar solution for r << tbh and breaks down before the PBH forms. Thus one still 
needs to consider both the initial perturbation shape and handle the formation of shocks. 
The contribution of these effects to the (initial state) entropy is likely to be small since the 
radiation density (and thus the entropy density) falls off at larger radii where self-similarity 
breaks down. 

For similar reasons including curved space effects will not induce an entropy much larger 
than computed here. For the metric 

ds^ = -a(r, tfde + A(r, tfdr^ + r^dn\ (Al) 

our expression for the "mass" of the radation (physical energy density integrated over the 
PBH volume) is 

Mrad = ^T^ p{r)A{r,tydr. (A2) 

This value is finite when computed for times before event horizon formation. At an infinites- 
imal amount of time before PBH formation, M^ad = ^pbh- The expression for the entropy 
is 

Srad = 47r X - i^^j j^ p{Tf''A{r, tydr. (A3) 

Compare these expressions to Equations ll|2|10l above. Because the entropy scales with den- 
sity as a power less than unity, non-uniform density profiles that maintain the normalization 



of Equation IA2I will tend to decrease Srad- 
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